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Abstract. Spectra of extranuclear clouds within the ionization 
cone of Circinus galaxy are presented and modelled assuming 
photoionization by the Seyfert 2 nucleus but no preconceived 
assumption concerning the spectral shape of the ionizing AGN 
radiation or the gas density distribution. 

The most important result is that, regardless of the assumed 
AGN spectral shape or density distribution, the metal abun- 
dances are remarkably well (±0.2 dex) constrained by the data. 
Oxygen and neon are found to be a factor of 5 below solar, 
while nitrogen is solar (cf. Table ||). The large N/O overabun- 
dance is in agreement with the chemical enrichment expected 
from the old (circum)nuclear starburst. Similar studies of ex- 
tranuclear clouds in other objects could therefore provide a 
powerful tool to determine metallicities and trace past starburst 
activity in AGNs. 

Other implications such as the role of shock excitation, 
which is effectively excluded, and the intrinsic shape of the 
AGN ionizing continuum, which is poorly constrained by the 
data, are also discussed. 
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1. Introduction 

The Circinus galaxy (A1409-65) is a nearby (~4 Mpc) gas 
rich spiral lying close to the galactic plane in a region of rela- 
tively low (Ay— 1.5 mag) interstellar extinction (Freeman et al. 



emission in polarized light (Oliva et al. 1998). Complementary 



1977). Several observed characteristics indicate that this galaxy 



hosts the nearest Seyfert 2 nucleus known. These include op- 
tical images sh owing a spectacular [OIII] cone (Marconi et al. 
1994 hereafter M94); optic al/IR spectra rich in prominent and 
narrow coronal lines (Oliva et al. 1994, hereafter 094 Moor- 
wood et al. 1996, hereafter M96); X-ray spectra displaying a 
very prominent Fe-K fluorescent line (Matt et al. 1996) and op- 
tical spectropolarimetric data which reveal relatively broad Ha 
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* Based on observations collected at the European Southern Obser- 
vatory, La Silla, Chile 



to these is observational evidence that this galaxy has recently 
experienced a powerful nuclear starburst which is now traced 



by the near IR emission of red supergiants (Oliva et al. 1995 



Maiolino et al. 1998) and which may have propagated outwards 



igniting the bright ring of O stars and HII regions visible in 



the Ha image (M94). Such a starburst could have been trig- 



gered by gas moving toward the nuclear region and eventually 
falling onto the accretion disk around the black hole powering 
the AGN. 

A debated issue is whether nuclear starbursts are common 
features of AGNs and if they are more common in type 2 
than in type 1 Seyfer ts, as suggested by e.g. 10/iin observa- 
tions (Maiolino et al. 1995) and studies of the stellar mass to 



light ratios ( 094 ). Since starbursts are predicted and observed 
to deeply modify the chemi cal ab undances of the host galaxy 
(e.g. Matteucci & Padovani 1993 ), such an effect should also 
be evident in this and other Seyferts. However, to the best of 
our knowledge, no reliable measurement of metallicity for the 
narrow line region clouds of Seyfert 2's exists in the litera- 
ture. In particular, although has been since long known that the 
large [NII]/Ha ratio typical of Seyferts cannot be easily ex- 
plained using simple models with normal nitrogen abundances 
(e.g. Osterbrock 1989} , Komossa & Schulz 1997), the ques- 
tion of whether its absolute (N/H) or relative (e.g. N/O) abun- 
dance is truly different than solar is still open. Finding a reli- 
able method to derive metallicities and, therefore, to trace and 
put constraints on past starburst activity is the main aim of this 
paper. 

We chose the Circinus galaxy as a benchmark because its 
emission line spectrum is characterized by remarkably narrow 



(Si 150 km/s, 094 ) emission lines which are particularly easy to 



measure and which indicate relatively low dynamical activity. 
This last aspect may be used to put tight constrains on the pos- 
sible contribution of shock excitation which may complicate 
the modelling of the observed spectrum and the determination 
of metallicities. 

This paper presents new optical and infrared spectroscopic 
data and is structured as follows. Observations and data reduc- 
tion are described in Sect. 2 and the results are analyzed in 
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Table 1. Observed and dereddened line fluxes 
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Fig. 1. Slit positions overlaid on grey-scale line images. The 
positions of the various knots are marked. 






Fig. 2. Spectrum of the nucleus, i.e. the central region at 
PA=318° (cf. Figs. 1, 4). Fluxes are in units of 10~ 16 erg cm~ 2 



erg cm s 



Sect. 3. In Sect. 4 we constrain the excitation conditions of the 
gas and model the observed spectra in terms of photoionization 
from the AGN. The derived chemical abundances are discussed 
in Sect. 5 and in Sect. 6 we draw our conclusions. 



2. Observations and Data Reduction 

Long slit optical spectra were collected at the ESO NTT tele- 
scope in March 1995 using a dichroic beam splitter feeding 
both blue (1024 2 Tek-CCD with 0.37'7pix) and red (2048 2 
Tek-CCD with 0.27'7pix) arms of EMMI. Simultaneous blue 
(3700-5000 with -1.7 A/pix) and red (5000-10000 with ~1.2 
A/pix) spectra were obtained through a 2"slit centered on the 
optical nucleus at two position angles (cf. Fig. [j]). The first was 
at PA=318°, roughly aligned with the [OIII] cone axis and in- 
cluding the brightest [OIII] emitting knots, while the second 
was at PA=243° and along the low excitation rim visible in the 
[SII] and "true colour" images shown in Fig. 5 and Fig. 10 of 
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Fig. 3. Spectrum of Knot C, i.e. a region 15.5"from the nucleus 
at PA=318° (cf. Figs. 1, 4). Flux units are as in Fig. 



M94j Q- Several exposures were averaged and the total inte- 
gration times of the spectra at PA=318° were 100, 75, 25 min- 
utes over the 3700-5000, 5000-7300 and 7700-10080 A wave- 
length ranges, respectively. The spectra at PA=243° covered 
3700-5000 and 5000-7300 A with a total integration time of 
50 minutes per wavelength interval. Data were flux calibrated 
using observations of LTT3218 and reduced with MIDAS us- 
ing standard procedures. 

Infrared spectra were collected in March 1995, also at the 
ESO NTT using the spectrometer IRSPEC equipped with a 
62x58 SBRC InSb array whose pixel size was 2.2"along the slit 
and ~5 A along the dispersion direction. Spectra of [Fell] 1 .644 
lira, H 2 2.121 /im, Br7 and [SiVI] 1.962 were obtained us- 
ing a 4.4"slit at PA=318°. Each long-slit spectrum consisted of 
4 ABBA cycles (A=source, B='sky', i.e. a region 300"E) with 
2x60 sec integrations per position. The data were reduced using 
the IRSPEC context of MIDAS which was developed by one of 
us (EO) and which also allows the accurate subtraction of time 
variable OH sky lines. The same airglow emission was used as 
wavelength reference (Oliva & Origlia 1992 ) and the spectra 
were corrected for instrumental and atmospheric transmission 
using spectra of featureless early O stars and flux calibrated 
using measurements of photometric standard stars. 

3. Results 

3.1. Line fluxes and ratios 

The quasi-complete spectra of the nucleus and of knot C are 
shown in Figs || respectively, and the derived line fluxes are 
summarized in Table [j]. Dilution by a stellar continuum is par- 
ticularly strong in the nuclear spectrum where the equivalent 



width of [OIII]A5007 is only 50 A (cf. Table ||) and a factor of 
~10 lower than found in typical Seyfert 2's. The stellar con- 
tribution is normally estimated and subtracted using either off- 
nuclear spectra extracted from the same 2D long slit frames, or 
a suitable combination of spectra of non-active galaxies used 
as templates (e.g. Ho 1996} , Ho et al. 1997). However, neither 
of the methods proved particularly useful because line emis- 
sion contaminates the stellar emission all along the slit, and 
we could not find any template which accurately reproduces 
the prominent stellar absorption features typical of quite young 
stellar populations. The fluxes of weak lines (<5% of the con- 
tinuum) in the nucleus are therefore uncertain and, in a few 



cases, quite different than those reported in 094 , the largest 
discrepancy being for [NI] which is a factor of 2 fainter here. 

The spectrum of knot C has a much more favourable line/- 
continuum ratio and shows many faint lines which are particu- 
larly useful for the modelling described in Sect. |j. 

3.2. Spatial distribution of emission lines 

The spatial variation of the most important lines is visualized 
in Figs. Q U which show contour plots of the continuum sub- 
tracted long slit spectra and selected spectral sections of the 
various knots respectively. The fluxes are summarized in Ta- 
ble U together with the extinctions which were derived from 
hydrogen recombination lines assuming standard case-B ratios 



(Hummer & Storey 1987). A remarkable result is the large vari- 



ations of the typical line diagnostic ratios [OBTJ/H/3, [OI]/Ha, 
[NII]/Ha and [SII]/Ha which are plotted in Fig. ^ and range 
from values typical of high excitation Seyferts (nucleus, knots 
A, B, C, D), to low excitation LINERs (knots H, I) and nor- 
mal HII regions (knots E, L). Another interesting result is the 
steep extinction gradient between the regions outside (knots C, 
D, H, I) and those close to the galactic disk (nucleus and knots 
A, E, L). How ever, a comparison betw een th e Br7 map (Moor- 
wood & Oliva 1994), the Ha images ( |M94) and the observed 
Bra flux from the whole galaxy ( |VI96 ), do not show evidence 
of more obscured {Ay ~ 10 — 30) ionized regions such as those 
observed in N GC494 5 and other starburst galaxies (e.g. Moor- 
wood & Oliva 1988). Nevertheless, these data cannot exclude 



the presence of deeply embedded ionized gas which is obscured 
even at 4/im (i.e. Ay>50 mag). 

Particularly interesting is the variation of the line ratios be- 
tween the adjacent knots C and D. The ratio [FeVII]/H/3 is a 
factor of 2 larger in knot D than in C, but this most probably 
reflects variations of the iron gas phase abundance (see also 
Sect. 4.5). Much more puzzling is the spatial variation of the 
low excitation lines [OI], [SII], [Nil] which drop by a factor 
1.8, while the high excitation lines Hell, [OIII], [Nelll], [Arlll] 
together with the [SII] density sensitive ratio and [OIII]/[OII] 
vary by much smaller amounts (cf. Table |J). This cannot there- 
fore be explained by variations of the ionization parameters, 
which should first of all affect the [OIII]/[OII] ratio. A possible 
explanation for this is discussed in Sect. p~8} 



These images are available at }ittp://www.arcetri.astro.it/^oliva 
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Table 2. Fluxes of significant lines in the various knots' 1 ) 
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• ' Dereddened fluxes relative to H/3=100, extinctions are computed imposing Hq=290 and the adopted visual extinctions are given at the top 

of each column. Blank entries are undetected lines with non-significant upper limits 
a H/3 is weak and the derived extinction is therefore uncertain. 

Possibly overestimated due to contamination by foreground gas with lower extinction. 
c Dereddened flux, units of 10 -15 erg cm -2 s -1 
d Equivalent widths in A of nebular emission lines. 



3.3. Diagnostic line ratios 

Temperature and density sensitive line ratios are summarized 
in Table g, Note that only a few values are available on the 
nucleus because of the strong stellar continuum which prevents 
the measurement of faint lines. 



3.3.1. Electron temperature 

The relatively large T e (OIII), slightly lower T e (SIII) and much 
cooler [Nil], [SII] temperatures are typical of gas photoion- 
ized by a "typical AGN", i.e. a spectrum characterized by a 
power law continuum with a super-imposed UV bump peaked 
at «50-100 eV (e.g. Mathews & Ferland |l987| ). As T e (OIII) 
mainly depends on the average energy of hv < 54 eV ion- 
izing photons, "bumpy" spectra, which are quite flat between 
13 and 54 eV, yield hot [OIII]. The lower ionization species, 
such as [Nil] and [SII], mostly form in the partially ionized re- 
gion, heated by soft X-rays, whose temperature cannot exceed 
10 4 K due to the powerful cooling by collisionally excited Lya 
and 2-photon emission. The contrast between the OIII and Nil 
temperatures can be further increased if sub-solar metallicities 
are adopted, because [OIII] is a major coolant while [Nil] only 
plays a secondary role in the cooling of the partially ionized 
region. 



An alternative explanation for the OIII, Nil etc. tempera- 
ture differences is to assume that part of the line emission arises 
from density bounded clouds. In this case there is no need to 
adopt a "bumpy" AGN spectrum and detailed models, assum- 
ing a pure power law ionizin g con tinuum, were developed by 
Binette et al. (1996, hereafter |B9(J). However, it should be kept 



in mind that T e (OIII)>r e (NII) does not necessarily indicate 
the presence of density bounded clouds. 

3.3.2. Electron density 

There is a clear trend between n e and excitation of the species 
used to determine the density. 

The [SII] red doublet yields densities lower than [Oil] 
which is compatible with a single-density cloud for the follow- 
ing reason. If the flux of soft X-rays (200-500 eV) is strong 
enough, then [SII] lines are mostly produced in the X-ray 
heated region where the average hydrogen ionization fraction is 
quite low (^0.1). The lines of [Oil] on the contrary can only be 
produced in the transition region, where the ionization degree is 
close to unity, because of the very rapid O-H charge exchange 
reactions. Hence n e (SII)<?i e (OII) most probably indicates the 
presence of a strong soft-X flux, as one indeed expects to be 
the case for an AGN spectrum. This is also confirmed by the 
detailed modelling described below. 

The higher densities in the fully ionized region are new re- 
sults because the blue [ArlV] doublet is usually too weak in 
AGN spectra and the FIR [ NeV] l ines are only accessible with 
the ISO-SWS spectrometer ( |M96| ). The [ArlV] density of knot 
C is equal, within the errors, to that derived by [Oil] thus in- 
dicating that no large variations of densities are present in this 
cloud. 



4. Photoionization models 

Several photoionization models for the Circinus galaxy have 
been discussed in the literature. These used line intensities and 
ratios measured from the central regions of the galaxy and were 
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Fig. 4. Intensity contour plots in the position-A plane. The long 
slit spectra are continuum subtracted and the levels are loga- 
rithmically spaced by 0.2 dex. The ordinate are arc-sec from 
the Ha peak along the two slit orientations (cf. Fig. [l]). The 
dashed lines show the regions where the spectra displayed in 
Figs. ||, |3|, |5] were extracted. 

Table 3. Temperature and density sensitive line ratios. 
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Value 
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Fig. 5. Spectra of the various knots (cf. Figs. |lj, |4|) at selected 
wavelength ranges including [OIII], H/3, Hell (left panels) and 
[FeVII], [OI], [Nil], Ha, [SII] (right hand panels). Fluxes are 
in units of 10~ 16 erg cm~ 2 s _1 A -1 and A's are in A. The 
spectra are also scaled by a factor given in the plots to show 
faint features. 



mostly directed to constraining the intrinsic shape of AGN pho- 
toionizing continuum. Quite different conclusions were drawn 
by M96, who found evidence for a prominent UV bump cen- 
tered at ~70 eV, and Binette et al. (1997, hereafter |B97b who 



showed that equally good (or bad) results could be obtained 
using a combination of radiation and density bounded clouds 
photoionized by a pure power-law spectrum. In all cases the 
metal abundances were assumed to be solar and no attempt was 
made to constrain metallicities. 

Here we mainly concentrate on knot C, an extranuclear 
cloud whose rich emission line spectrum and simple geometry 
(a plane parallel region) are better suited for deriving physical 
parameters from photoionization modelling. We first analyze 
the observational evidence against shock excitation and then 
describe in some detail the new modelling procedure which is 
primarily aimed at determining the gas metal abundances, (cf. 



Sects. 4.2 and 5.1). We also analyze the problems to repro duce 



the observed [FeVII]/[FeII] and [NII]/[NI] ratios (Sects. 4.5 



]), discuss the role of dust and reanalyze several crucial as- 
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Fig. 6. S pectr al classification diag rams (from Veilleux & Os- 
terbrock 1987 and Baldwin et al. 1981 ) using the line ratios 
observed in the various knots. Note that "N" refers to the nu- 
clear spectrum while "I" is the average of the H and I knots 
which have similar line ratios (cf. Table 2). 



pects of the line emission from the brightest regions closest to 
the nucleus. 



4.1. Arguments against shock excitation 

Strong evidence against shock excitation comes from the very 
low strength of [Fell] which indicates that iron is strongly un- 
derabundant, i.e. hidden in grains, within the partially ionized 
region (cf. Figs. |l0[ [H]). Shocks are very effective in destroy- 
ing dust grains, and velocities as low as 100 km/s are enough to 
return most of the iron to the gas phase (e.g. Draine & McKee 
This is confirmed by observations of nearby Herbig- 



1993) 



Haro objects (e.g. Beck-Winchatz et al. 1994) and supernova 



remnants whose near IR spectra are dominated by prominent 



[Fell] IR lines (e.g. Oliva et al. |19890 . Although the low Fe+ 
gas phase abundance could, in principle, be compatible with a 
very slow shock (^50 km/s), this falls short by orders of mag- 
nitude in producing high excitation species. 

Another argument comes from the HeII/H/3 and Hell/Hel 
ratios which are a factor >2 larger than those predicted by 
shocks models with velocities v < 500 km/s (Dopita & Suther- 
land 



1995) 



It should be noted that, although stronger Hell 
could probably be obtained by increasing the shock speed to 
^1000 km/s, these velocities are incompatible with the ob- 
served line profiles. More generally, the observed line widths 



(<150 km/s, 094) are difficult to reconcile with the fact that 
only shocks faster than 300 km/s can produce promine nt hig h 
excitation lines such as [OIII] (Dopita & Sutherland 1996). 



This argument becomes even stronger when interpreting the 
nuclear spectra where the highest excitation lines, [FeX,XI], 
do not show any evidence of large scale motions at velocities 
>200 km/s. (cf. p94 and work in preparation). 

Finally, it is worth mentioning that detailed shock + pho- 
toioni zation composite models recently developed by Contini 
et al. ( 1998 ) suggest that shock excited gas does not contribute 
significantly to any of the optical/IR lines from Circinus. 

4.2. Details of the photoionization models 

4.2. 1 . Single component, radiation bounded clouds 

We constructed a grid of models covering a wide range of ion- 
ization parameters (—3.5 < log U < —1.5), densities (2.0 < 
log n < 5.0), metallicities (0.08 < Hc/H < 0.16, -1.5 < 
log Z/Zq < 0) and shape of the ionizing continuum which we 
parameterized as a combination of a power law extending from 
1 eV to 10 keV with index 2 < a < 0.4 and a black-bodyf] 
with 5 < XogTsB < 6, the relative fraction of the two com- 
ponents being another free parameter. All the parameters were 
varied randomly and about 27,000 models were produced using 



Cloudy (Ferland 1993) which we slightly modified to include 
the possibility of varying the N-H charge exchange rate. The 
line intensities were also computed using the temperature and 
ionization structure from Cloudy and an in-house data base of 
atomic parameters. The agreement with the Cloudy output was 
good for all the "standard" species while large discrepancies 
were only found for coronal species (e.g. [FeX]) whose colli- 
sion strengths are still very uncertain and much debated. 

Out of this large grid we selected about 400 models 
whose [OIII]/H/3, [OIII]/[OII]/[OI], [ArV]/[ArIV]/[ArIII], 
[SIII]/[SII], [OIII]A5007/A4363, [OII]AA3727/AA7325, 
[SII]A6731/A6716, [ArIV]A4711/A4740 line ratios were 
reasonably close to those observed in knot C. These were used 
as the starting point for computing the "good models", with 
adjusted values of relative metal abundances, which minimized 
the differences between predicted and observed line ratios. 
Note that to reproduce both [FeVII] and [Fell] we were in all 
cases forced to vary the iron gas phase abundance between 
the fully and partially ionized regions (cf. also Sect. _p~5| ). 
The results of the best model are summarized in Table |j and 
discussed in Sect. 43. Note the large discrepancy for [NI] 
which is overpre dicte d by a factor of ^5, this problem is 



discussed in Sect. 4.6 



The most important results are the abundance histograms 
shown in Figs. [l(], [jl] which were constructed by including all 
the "good" models with Xred <• Although the choice of the 
cutoff is arbitrary, it should be noticed that variations of this 
parameter do not alter the mean values, but only influence the 



The choice of a black-body does not have any direct physical 
implications. We used it just because it provides a simple way to pa- 
rameterize the UV bump which seem s to be a typical features of AGNs 
(e.g. Mathews & Ferland 



1987 



Laor 

tion disk models (e.g. Laor & Netzer 



1989) 



1990 ) and is predicted by accre- 
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shape of the distributions whose widths roughly double if the 

^rea 



2 , cutoff is increased to values as large as 30. 



4.2.2. Multi-density components, radiation bounded clouds 

The large grid described above was also organized to have at 
least 4 models with the same photon flux, continuum shape 
and abundances but different densities, and these were used 
as starting points to construct multi-density models. We sim- 
ulated 2-density clouds by coupling all available models with 
different weights and selected about 300 of them, which were 
used as starting points to compute the "good models" follow- 
ing the same procedure adopted for the single-density case. We 
also simulated clouds with >3 density components, but these 
results are not included here because this complication had vir- 
tually no effect on the quality of the fit. 

The most important result is that single and multi-density 
models are equally good (or bad) in reproducing the observed 
line ratios. Obviously, this does not necessarily imply that knot 
C is a single-density cloud, but rather indicates that the stratifi- 
cations which probably exist have little effect on the available 
density sensitive line ratios. 

4.2.3. Mixed models with density and radiation bounded 
clouds 

We constructed a grid of about 10,000 models photoionized 
by a power law AGN continuum with index 2 < a < 0.3, which 
turns into j/~ - 7 beyond 500 eV, and covering the same range of 
physical parameters as the radiation bounded clouds described 
above. The column density of the radiation bounded compo- 
nent was always large enough to reach a temperature <3000 
K at its outer edge. For each model we also computed line in- 
tensities from a density bounded cloud which we defined as a 
layer with thickness A R equal to 5% of the nominal radius of 
the HII Stromgren sphere, numerically (/= filling factor) 



Ai?~ 2500 



U 
nn f 



pc 



The relative contribution of density and radiation bounded re- 
gions was set by imposing HeIIA4686/H/3=0.6 



This approach is somewhat similar to that adopted by B96 



apart from the following details. The radiation bounded compo- 
nent here is formed by high density gas, a condition required to 
match the observed [ArIV]A471 1/A4740, and is concentrated 
within the projected size of knot C, i.e. <80 pc. In practice, the 
two components have the same densities and see un-filtered 
radiation from the AGN. 

The "good models" were optimized, and the best abun- 
dances derived using the the same procedure depicted above. 
Worth mentioning is that most of the "good models" have 
spectral slopes in the range 1.3-1.5 which are in good agree- 



ment with those used by B96 The most important result is that 



4.2.4. The role of dust in extranuclear knots 

Dust can modify the ionization and temperature structure be- 
cause it competes with gas in absorbing the UV photons, and 
because it hides refractory elements such as iron. Given the low 
ionization parameters, however, the first effect is negligible, i.e. 
the ionization structure of the fully ionized region of knot C is 
not affected by dust although this plays an important role in 
modifying the heating-cooling balance of the partially ionized 
region heated by X-rays. 

Cooling: the refractory species hidden in the grains cannot 
contribute to the line cooling and this effect is correctly com- 
puted by Cloudy. For example, depleting Fe on grains produces 
higher gas temperature and stronger [OI], [SII] etc. lines, be- 
cause it suppresses the near IR lines of [Fell] which are among 
the major coolants for quasi-solar Fe/O gas phase abundances. 

Heating: the metals hidden in the dust still contribute to 
the heating of the gas because the X-rays have energies much 
larger than the binding energy of the grains, and cannot there- 
fore recognize metals in dust from those in the gas phase. 

The major problem is that Cloudy (and probably other pho- 
toionization models) does not include the X-ray heating from 
the grain metals and therefore underestimates the temperature 
of the partially ionized region and hence may predict weaker 
fluxes of [OI], [SII] and other low excitation lines. Therefore, 
the models for knot C, having most of iron depleted on grains, 
are not fully self-consistent and most probably require a too 
high flux of X-rays to reproduce e.g. [OI]/[OHI]. This may im- 
ply that the "true" models should have somewhat softer spectra 
than those of Fig. ^. 

4.3. Best photoionization models of knot C: the shape of the 
AGN continuum and the role of density bounded clouds 

Table lists the results of the best photoionization models, se- 
lected from the several hundred which provided a reasonable fit 
to the spectru m of k not C. Note that the results of multi-density 
clouds (Sect. [4.2.2 ) are not included because they are virtually 
indistinguishable from those of single-density models. 

Fig. ^ shows the AGN continuum adopted for the "best 
models" of Table [I| Radiation bounded models require a 
prominent UV bump centered (in F y ) at about 100 eV. This 
is slightly bluer than that found by M96 and somewhat harder 



than model spectra of more luminous accretion disks (cf. e.g. 



Laor 1990), but in qualitative agreement with the predicted 



dependence of spectral shape with AGN luminosity, i.e. that 
lower luminosity accretion disks should have harder spectra 



(Netzeretal. 1992) 



these mix ed models give similar, though slightly poorer (cf. 
Sect. 4.3 ) results than radiation bounded clouds. 



The amplitude of the UV bump could be significantly 
decreased by relaxing the assumption on the gas density 
distribution and, with a properly tuned combination of 
density and radiation bounded clouds, one may probably 
get a similarly good fit with a power law. It is nevertheless 
instructive to analyze why our mixed models provide a worse 
fit to the data and, in particular, are unable to simultaneously 
reproduce the OII/OIII balance and the [ArV]/[ArIII] and 
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Table 4. Photoionization models for knot C 



Adopted parameters 


Model 1 


Model 2 


Model 3 


Element 




Abundances' 1 - 1 




Helium 




11.1 


11.1 


11.1 


Nitrogen 




8.05 


8.20 


8.05 


Oxygen 




8.17 


8.27 


8.17 


Neon 




7.27 


7.37 


7.27 


Sulphur 




6.98 


6.90 


6.75 


Argon 




6.38 


6.45 


6.20 


Silicon 




6.85 


6.95 


U.OJ 


Fe low-ion a 




5.80 


5.80 


5.80 


Fe high-ion 6 




6.90 


7.00 


7.00 


Gas density (uh) 


2000 


3000 


ZUUU 


Ionizing continuum c 


UV-bump 




V 


Ionization parameter d (U) 


3.5 10~ 3 


4.7 10~ 3 




Fraction radiation bounded 6 


100% 


4% 


3% 


F(KP)/F RB (H/3) f 


100% 


6% 


6% 


Line ratio' 2 ) 


Observed 


MnHpl 1 
1V1UUC1 1 


MnHpl ? 


Model 3 


HeII/H/3 


0.60±0.08 


U.oO 


O.oU 


0.60 


Hell/Hel 


6.2±0.9 


u.o 


J.O 


6.3 


[NII]/Ha 


0.80±0.08 


0.83 


0.79 


0.77 


[NII]/[NI] 


34±7 


8.0(!!) 


6.5(!!) 


4.7(H) 


A6583/A5755 


50±13 


41 


57 


39 


[OIII]/H/3 


9.7±1 


10 


9.8 


9.3 


A5007A4363 


46±9 


31(!) 


43 


42 


[OII]/[OIII] 


0.14±0.02 


0.12 


.057(1!) 


0.12 


[OI]/[OIII] 


.027±.003 


.027 


.024 


.028 


AA3727/AA7325 14±4 




1 z. 


13 


[NeIII]/H/3 


0.66±0.1 


0.68 


0.70 


0.73 


[SiVI]/Br 7 


<5 


2.2 


1.5 


0.12 


[SII]/Ha 


0.19±0.02 


0.20 


0.20 


0.19 


A6731/A6716 


0.89±0.06 


0.95 


0.88 


0.91 


[SIII]/[SII] 


2.8±0.4 


3.0 


2.4 


2.6 


A9531/A6312 


26±7 


20 


22 


21 


[ArIII]/Ha 


.092±.011 


.089 


.098 


.092 


[ArIV]/[ArIII] 


0.37±0.1 


0.56 


0.57 


0.25 


[ArV]/[ArIII] 


0.24±0.03 


0.17 


0.17 


.032(H) 


A4711/A4740 


1.0±0.2 


0.88 


0.96 


0.88 


[FeVII]/Ha 


.033±.005 


.031 


.030 


.005(H) 


[FeVI]/[FeVII] 


<0.7 


0.70 


1.2 


4.1(H) 


[FeX]/[FeVII] 


<0.2 


.045 


.030 


.005 


[FeII]/Br7 


wl.5 


1.6 


1.8 


1.8 



^ 12+log(X/H) where X/H is the absolute abundance by number 
(2) Abbreviation used: HeII=A4686, HeI=A5876, [NII]=A6583, 
[NI]=AA5200, [OIII]=A5007, [OII]=AA3727, [OI]=A6300, 
[NeIII]=A3869, [NeV]=A3426, [SiVI]=A 19629, [SiVII]=A24827, 
[SII]=A6731, [SIII]=A9531, [ArIII]=A7136, [ArIV]=A4740, 
[ArV]=A7006, [FeVII]=A6087, [FeVI]=A5146, [FeII]=A16435 

a Gas phase abundance in the partially ionized region 

b Gas phase abundance in the fully ionized region 

c UV-bump has T BB =4 10 5 K and a=1.4, cf. Fig. | 

d U = Q(H)/4TYR 2 n H c 

e Models 2,3 include density bounded clouds with the same density 

and U as radiation bounded regions, cf. Sect. 4.6.2 
' Flux relative to that produced by a radiation bounded model 
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Fig. 7. The AGN ionizing continua used to compute the best 
models (Table |j) are plotted together with the observed X and 
IR continua. The solid line refers to a single density, radia- 
tion bounded component which provides a good fit for all lines 
but [NI]. The broken curves are for combinations of density 
and radiation bounded clouds, but these models cannot simul- 
taneously reproduce the [OII]/[OIII] and [ArV]/[ArIII] ratios. 
The dashed curve provides the best fit for high excitations lines 
while the dotted curve best reproduces the OI/OII/OIII ioniza- 
tion balance (see Sect. 4.3 for details). Note that the curves 



show the spectra seen by knot C, and should be scaled by a 
coefficient which takes into account the "intrinisic beaming" 
of the AGN ionizing radiation (a factor of 2 for an optically 
thick disk) before being compared with the observed points. 



See Sect. 4.4 for a discussion of the AGN energy budget. 



[FeVII]/[FeVI] ratios (cf. last two columns of Table Q). Model 
#2 has a high ionization parameter and correctly predicts 
high excitation lines but underestimates all [Oil] lines while 
model #3, with a lower value of U, correctly reproduces the 
[OII]/[OIII] ratio but predicts very faint high excitation lines. 
The main reason for the above difference is that, for a given 
gas density, the [OII]/[OIII] ratio is a measure of the flux of 
soft ionizing photons (hv= 13-54 eV), while [ArV]/[ArIII] and 
[FeVII]/[FeVI] depend on the flux of harder radiation {hv <; 6 
80 eV). The spectrum of knot C is characterized by a quite 
large [OII]/[OIII] ratio, which points towards low fluxes of 
soft photons, and strong high excitation lines, which require 
large fluxes of hard photons. These somewhat contradictory 
requirements can be easily satisfied by a spectrum which 
steeply rises beyond the Lyman edge, and peaks at «100 eV, 
i.e. a spectrum similar to that of Model #1. 

Independently on the detailed results of the models, the fol- 
lowing arguments indicate that density bounded clouds may in- 
deed play an important role. 

- Ionizing spectra with pronounced UV-bumps tend to pro- 
duce too high [OIII] temperatures and model #1 is indeed quite 
hot, though still compatible (within 1.5er) with the somewhat 
noisy measurement of [OIII]A4363 (cf. Fig. || and Tab. Q). Cu- 
riously, though, the result of Model #1 is the opposite of that 
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Table 5. Global properties of the AGN ionizing continuum^ 1 ) finds Q{9) w cos 9 (e.g. Laor & Netzer 1989| ) and 



Ionizing flux seen by knot C a (cm 2 s 1 ) (1-4) x 10 11 

Effective covering factor of knot C b ~5% 
Q(H)i GN ' (s" 1 ) (0.5-2.0) x 10 54 

Observed recombination rate d 
Fraction of Q(H)agn intercepted by gas 



(s" 1 ) 2 ■ 10 52 

~1% 



Total AGN ionizing luminosity-' 
Observed FIR luminosity (Lfir) 9 



(Lq) 



~ 2 ■ 10 1 
1.2 ■ 10 1 



^ Assuming a distance of 4 Mpc 

a Required to reproduce the I /-sensitive a n d d e nsity sensitive line 

ratios measured in Knot C (cf. Sects. |4.3| , 
b Ratio between the covering factor and the ~ 40 x 40 pc 2 geometric 

cross section of the knot 
c Total number rate of AGN ionizing photons, assuming emission 

from an optically thick disk (cf. Sect. }AM 
d From L(Bra)=2.5 10 5 L© (M96) 

e Assuming that half of the observed Bra emission is produced by 

the starburst ring 
' Assuming an average phot on ene rgy of 50 eV 
9 From Siebenmorgen et al. (1997) 



obtained by most previous models in the literature which failed 
to produce hot enough [OIII]. Therefore, the classical "[OIII] 
temperature problem" may just reflect the fact that past mod- 
els were mostly biased toward high metallicities and rather fiat 
(i.e. without strong bumps) continua. 

- The H/3 flux from knot C is only ~5% of that expected 
if the gas absorbed all the ionizing photon impinging on the 
2" x 2"(40 x 40 pc 2 ) geometrical cross-section of the cloud. 
Such a small "effective covering factor" could, in principle, 
be obtained by assuming a suitable distribution of radiation 
bounded, geometrically thin clouds or filaments. However, den- 
sity bounded clouds seem to provide a more self-consistent in- 
terpretation because the measured value is remarkably close to 
the 6% predicted by models #2 and #3 (cf. note / of Tab. Q). 

- The variation of line ratios between the different knots can- 
not be explained by radiation bounded clouds illuminated by 
the same ionizing continuum, but require e.g. intrinsic beam- 
ing of the AGN continuum and/or filtering by matter bounded 
clouds closer to the nucleus. Mixed models could give a more 
natu ral e xplanation to the spatial variation of line ratios (see 
also 5g5). 



4.4. Energy budget of the AGN 

The ionizing photon flux seen by knot C is constrained by the 
[/-sensitive and density sensitive line ratios which yield 



U-nu -c~ (1 - 4) x 10 11 cm 



This can be translated into Q(H)agn, the total number rate of 
ionizing photons from the AGN, once the angular distribution 
of the UV ionizing radiation is known. Assuming that it arises 
from a geometrically thin, optically thick accretion disk, one 



0(H) 



AGN 



2irR 



knC ^ion 



(0.5 - 2.0) X 10 54 _ x 



where i?knC is the projected distance of knot C i.e. 15.5"or 
e?=300/cos i pc from the nucleus, i being the inclination angle 
of the cone relative to the plane of the sky. Detai led kinematical 
studies indicate \i\ < 40° (Elmouttie et al. Il998|). 



The AGN luminosity in the ionizing continuum is therefore 



Li, 



Q(H) 



AGN 



< hv ion > 



1-4-10 



10 



-L, 



9 ■ ^0 

cos^ I 



where < hvi on > is the average photon energy which is here as- 
sumed to be ~50 eV. The ionizing luminosity is therefore very 
large but compatible, within the errors, with the observ ed FIR 
luminosity Lfir — 1.2 10 10 i© (Siebenmorgen et al. 1997 ). 
This implies that the AGN emits most of its energy in the ioniz- 
ing continuum or, equivalently, that the AGN intrinsic spectrum 
has a prominent peak or bump in the ionizing UV, and much 
weaker emission at lower energies. This is in good agreement 
with computed models of accretion disks which also predict 
that low luminosity AGNs, such as Circinus, should be char- 
acterized by a quite hard ionizing continuum (cf. Laor 1990 
Netzer et al 



1992) 



The global properties of the AGN are summarized in Ta- 
ble H which also includes a comparison between Q(H)agn and 
the observed recombination rate, based on IS O obs ervations of 
the Bra H-recombination line at 4.05 fxm (cf. |M96| ). The differ- 
ence is remarkably large, with only ~1% of the AGN ionizing 
photons being accounted for by emission from "normal" ion- 
ized gas. This indicates that the bulk of the Lyman continuum 
radiation either goes into ionizing regions which are obscured 
even at 4 /im (i.e. Av>50 mag), or is directly absorbed by 
grains in dusty clouds lying very close to the AGN. 

4.5. Iron depletion and the [FeVII]/[FeII] problem 

The observed [FeVII]A6087/[FeII]A16435^1 ratio cannot be 
explained using the same iron gas phase abundance in the Helll 
Stromgren sphere, where FeVII is formed, and in the partially 
ionized region, where iron is predominantly Fell due to the 
rapid charge exchange recombination reactions with neutral 
hydrogen. It should be noted that this is a fundamental prob- 
lem unrelated to the details of the photoionization models and 
primarily reflects the fact that the [FeVII] line has a very small 
collision strength (T/cJi ~ 0.1) which is factor of about 10 
lower than the [Fell] transition. A Fe +6 /Fe + il 2 integrated rel- 
ative abundance is thus required to reproduce the observed line 
ratio. This number is incompatible with the relative sizes of the 
Helll and partially ionized regions which are constrained by 
e.g. the Hell and [OI] lines. It is also interesting to note that this 
problem is even exacerbated if the shock models of Dopita & 



Sutherland (1996) are adopted because these predict Fe +6 /Fe H 



ratios much smaller than the photoionization models described 
above. 
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Fig. 8. Effect of varying S(N), the rate coefficient of N-H 
charge exchange, on the [NII]/[NI] and [NI]/H/3 ratios pre- 
dicted by the photionization models discussed in Sect. \3. The 
solid line refers to model #1 while the dashed curve is for 
model #2 (cf. Table f§), and the dashed area show the value 
(with errors) measured in knot C. The parameter <5o(N) is the 
standard rate coefficient used by Cloudy. 



A possible solution could be to advocate that the [FeVII] 
collision strengths are underestimated by a factor of ~10 which 
would reconcile the observed [FeVII] and [Fell] intensities 
with a low (2% of solar) iron abundance (cf. Figs. [h], [H]). How- 
ever, although the collision strengths of coronal Fe lines are 
known to be very uncertain, and vary by factors >10 depend- 
ing on whether resonances are included in the computation (cf. 
Sect. 5 of Oliva 1997), the available collision strengths for [Fe- 



VII] increase by only 20% between the old DW computations 



of Nussbaumer & Storey (1982) and the newer R-MAT (i.e. 



including resonances) values of Keenan & Norrington (1991). 
More detailed modelling of the [FeVII] lines and spectroscopic 
studies of nearby astrophysical laboratories (e.g. high excita- 
tion planetary nebulae) are required to clarify this issue. 

The alternative possibility is to assume that the iron deple- 
tion is much larger in the partially than in the fully ionized 



region, as already suggested by Netzer & Laor (1993) and Fer 



guson et al. ( 1997). However, we could not find any simple ex- 
planation for such a stratification in knot C which lies far from 
the AGN and receives a relatively modest flux of hard UV pho- 
tons. Therefore, Fe-bearing dust cannot be photo-evaporated 
and the only mechanism to destroy grains is sputtering. The 
slow shock produced by the ionization front which propagated 
outward when the AGN turned on was too slow (<40 km/s) 
to e ffectiv ely return Fe to the gas phase (e.g. Draine & Mc- 
Kee 1993| ). Faster shocks (>100 km/s) are a natural and effi- 



cient source of sputtering but cannot explain the observations 
because they are expected, and observed, to emit prominent 
[Fell] lines from the dust-free recombining gas. A possibil- 
ity to overcome this problem is a combination of shocks and 
photoionization where e.g. the dust-free gas processed by the 
shock is kept fully ionized by the AGN radiation. However, this 
situation is short lived because, after a few thousand years, the 
gas piled up behind the shock will eventually reach a column 
density high enough to become radiation bounded, and shield 
the recombining gas which will therefore show up in [Fell]. 

4.6. The [NII]/[NI] dilemma 

The photoionization models of knot C systematically overes- 
timate by large factors (>6) the strength of [NI] relative to 
[Nil] (cf. Table | and Fig. |l0|). Although [NI]AA5200 has 
a quite low critical density (~1500 cm -3 ), this error can- 
not be attributed to the presence of higher density clouds be- 
cause these would also effect the [SII] density sensitive ra- 
tio. In other words, multi-density models which correctly re- 
produce the high [NII]/[NI] ratio inevitably predict too large 
[SII]A6731/A6716 ratios. Also, we can exclude observational 
errors because, in knot C, the [NI] doublet has an equivalent 
width of about 2.5 A and is only marginally affected by blend- 
ing with neighbouruing stellar absorption lines. 

A possible solution is to argue that the rate coefficients for 
N-H charge exchange are overestimated, as already suggested 



by e.g. Ferland & Netzer (1983). In the partially ionized re 



gion, Nil is mostly neutralized via charge exchanges with H° 
and adopting lower charge exchange efficiencies yield larger 
[NII]/[NI] ratios. This is evident in Fig. ^] where this ratio is 
plotted as a function of the assumed value of S(N), the charge 
exchange rate coefficient. Assuming a N-H charge exchange a 
factor of ~30 lower than presently adopted yields the correct 
[NII]/[NI] ratio. 

Noticeably, the problem we find here is exactly the opposite 



of what reported by Stasinska (1984) whose models systemat- 



ically underpredicted the [NI]/[NII] ratio in objects with low 
[OI]/[OII]. 

4.7. Modelling the nuclear spectrum: dusty, dust-free and 
diffuse components. 

A puzzling aspect of the line emission from regions very 
close to the nucleus is the simultaneous presence of high (e.g. 
[FeXI]) and low (e.g. [SII], [OI]) ionization species. More 
specifically, the images of M94 clearly show that [SII] peaks 
at a distance of only ~0.5" or 10 pc from the nucleus (cf. Fig. 



9 of [VI94 ). This result is incompatible with the standard idea 
according to which low excitation lines are produced in clouds 
with low ionization parameters. Numerically, combining the 
ionizing continuum inferred from the spectrum of knot C with 
n=1.2 10 4 cm -3 , the highest density compatible with the FIR 
[NeV] doublet (Table yields U ~ 0.5 at r=10 pc from 
the AGN, an ionization parameter far too large for the produc- 
tion of low excitation species. Not surprisingly therefore, all 
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Fig. 9. Ionization structure and spatial variation of the line 
emission from two nuclear clouds exposed to the AGN contin- 
uum of "Model #1" (cf. Fig. ^ and with identical physical pa- 
rameters (riH=10 4 cm~ 3 and filling factor /=0.025) except for 
dust. The model without dust emits most of the coronal lines, 
while the dusty cloud accounts for the prominent low excita- 
tion lines observed close to the nucleus. However, both models 
fail to produce enough [OIV] and other intermediate ionization 
lines which therefore require a third, more diffuse component. 



See Sect. 4.7 for details. 



the models so far developed for the high excitation lines (094, 
M96, B97) predict that [SII] should peak much farther out and 



have a much lower surface brightness than that observed. 

A simple and indeed natural solution is to assume that the 
low excitation lines are formed in dusty clouds. At these large 
ionization parameters dust dominates the absorption of UV 
ionizing photons and, therefore, quenches the Helll and HII 
Stromgren spheres. Consequently, the X-ray dominated par- 
tially ionized region starts at much smaller radii, and is also 
slightly hotter than in the dust free model (cf. Fig ||). Therefore, 
dusty clouds have a [SII] peak much closer to the nucleus and 
a > 10 times higher surface brightness than dust-free clouds. It 
should be noted, however, that the total luminosity of low ex- 
citation optical lines is similar in the two cases, fundamentally 
because the available luminosity of soft X-rays is the same. 

Although the combined emission of dusty and dust-free 
clouds can account for the observed emission of low ionization 
and coronal species, it falls short by large factors in producing 
[OIV] and other relatively low ionization species which form 
within the Helll sphere. This is an intrinsic limitation of "com- 
pact models" such as those of Fig. ^, and can be understood as 
follows. Compact models are characterized by large "ionization 
parameters" (cf. Sect. 3 of Oliva 1997 for a critical discussion 



on this parameter) and therefore have large fluxes of OIV ion- 
izing photons which keep oxygen in higher ionization stages 
(mostly OVI and OVII) at all radii inside the Helll Stromgren 
sphere. Outside the Helll region, on the contrary, oxygen can- 
not be ionized above OIII because most of the OIII ionizing 
photons have already been absorbed by Hell. Therefore, OIV 
can only exist in a very narrow range of radii, just at the edge 



of the Helll sphere, and its relative abundance is therefore very 
low. 

In practice, we found it impossible to construct a single 
model which simultaneously produces a com pact coronal line 
region, such as that observed in [FeXI] ( 094 ), and which comes 
anywhere close to the [OIV]A25.9/[OIII]A5007~0.3 observed 
ratio. We did indeed construct many thousands of randomly 
generated dusty and dust-free models, and attempted an ap- 
proach similar to that used for knot C (Sect. 4.2 ) but, in no case, 
could we find a model which satisfies these contradictory con- 
straints. It should also be noted that B97 independently come 
to a similar conclusion. 

The main conclusion therefore is that, regardless of the de- 
tails of the models, the nuclear spectrum and line spatial distri- 
bution can only be modeled by adding a third "diffuse" com- 
ponent (e.g. with a lower filling factor) to the dusty and dust- 
free clouds discussed above and depicted in Fig. ^| Given the 
large number of free parameters, we abandoned the idea of us- 
ing photoionization models to constrain abundances and other 
physical properties of the gas, We made some attempt to verify 
that a mixture of clouds exposed to the same continuum, and 
with the same abundances as Model #1 of knot C (Table |] and 
Fig. ^) could reasonably well reproduce the observed proper- 
ties of the nuclear spectrum. However, the results are not too 
encouraging and, apart from the much improved [SII] surface 
brightness and [OIV]/[OIII] ratio, the solutions are not signif- 



icantly better than those already discussed by M96 and B97 
and are not therefore discussed here. 



4.8. Modelling other extra-nuclear knots 

An analysis similar to that used for knot C was also applied to 
the other extra-nuclear knots using the more limited number of 
lines available in their spectra. The results can be summarized 
as follows. 

The abundances derived in the Seyfert-type knots A, B, D, 
G, F (cf. Fig. g) are similar those found in knot C but affected 
by much larger errors because of the more limited numbers of 
lines available for the analysis. In particular, the density sensi- 
tive [ArlV] doublet and the U-sensitive [ArV] line is not de- 
tected in any of these knots, and the reddening correction for 
[Oil] AA3727 could be very uncertain in the high extinction re- 
gions (cf. note b of Table ||). 

We also attempted to verify if the observed line ratios in 
these knots could be explained as photoionization by the same 
AGN continuum seen by knot C but could not find any satis- 
factory solution using radiation bounded clouds exposed to the 
same continuum. Adding matter bounded clouds alleviates the 
problem (as already stressed by B96| ) but requires an ad hoc 
choice of their photoelectric opacities, i.e. the radius at which 
the ionization structure is cut. In particular, explaining the drop 
of low excitation lines between the adjacent knots C and D (cf. 



Sect. 3.2) requires matter bounded clouds cut at about 1.2x 



the HII Stromgren radius, the exact position of the cut depend- 
ing on the assumed shape of the AGN continuum. Another pa- 
rameter affecting the ratio of low-to-high excitation lines (e.g. 
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Fig. 10. Element abundances (in log of solar units) as derived 
from a comparison between the observed line strengths in knot 
C, and the predictions of a large grid of randomly generated 
photoionization models. The three columns refer to clouds with 
different gas density distributions, and the "good models" are 
those coming closer to reproducing the observed line ratios. 
See text, Sects. 5.1 and X.2 for details. 



[OI]/[OIII]) is the iron gas phase abundance which influences 



the cooling of the partially ionized region (cf. Sect. X.2 A). If 



iron is more abundant in knot D, as indicated by its stronger 
[FeVII] line emission (cf. Table and Sect. O), than [OI], 



[SII], [Nil] could be depressed by the increased [Fell] cooling. 

The abundances derived for the LINER-like knots (H, I) 
are very uncertain (±1.3 dex at least). Their spectra are not 
compatible with illumination from the same continuum seen 
by knot C but require a harder (i.e. more X rays relative to 13- 
80 eV photons) spectrum which could be in principle obtained 
by filtering the AGN continuum through an absorber with a 
carefully tuned photoelectric opacity. Alternatively, the weak 
(in surface brightness) spectrum of these low excitation knots 
could be explained by shock excitation, in which case one ex- 
pects [FeII]A16435/H/3~l and a factor >10 larger than in the 
case of pure photoionization. 

Finally, the oxygen lines in the highly reddened Hll-like 
knots (E, L) are too weak to allow any reliable abundance anal- 
ysis. 

5. Discussion 

5.7. Element abundances 

Deriving abundances of AGN clouds from photoionization 
models has generally been considered to be unreliable because 
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Fig. 11. Same as Fig. |l0| but using photoionization models 
where the rate coefficient for N-H charge exchange reactions 
were arbitrarily set to zero. Note that this provides a much bet- 
ter match for the [NII]/[NI] line ratio while has a small effect 
on the Nit rogen abundance derived from [Nil] lines. See text, 
Sects. 5.1 and X.2 for details. 



the shape of the AGN ionizing continuum and the density dis- 
tribution of the emitting regions (clouds) are both basically 
unknown. The more typical approach therefore has been to 
assume a metallicity and use photoionization models to con- 
strain the AGN spectrum and/or the gas density distribution, 
or just to demonstrate that the gas is photoionized. Although 
explicit statements concerning the nitrogen abundance are of- 
ten found in the literature (e.g. Storchi-Bergmann & Pastoriza 



1989, Simpson & Ward 1996 1 these are based on a very lim- 



ited choice of photoionization model parameters and in most 
cases find oxygen abundances close to solar, in disagreement 
with what is derived here. The only other piece of work which 
covers a model parameter range comparable to that presented 



here is that by Komossa & Schulz (1997) who analyzes a much 



more limited numbers of lines, e.g. do not include [ArIV,V], in 
a large sample of Seyferts. They find that, on average, oxygen 
is underabundant by a factor of ^2 and that the N/O ratio is 
only a factor of 1.5-2.0 above the solar value. We believe that 
the analysis presented here leads to more reliable metallicity 
estimates. 

The results presented here indicate that, in spite of the 
above uncertainties, reliable metallicities can indeed be derived 
from spectra including a large enough number of lines, such as 



that of knot C. Our method has been described in Sect. X.2 and, 



in short, is based on the computation of a large number of pho- 
toionization models with a minimal personal bias and precon- 
ceived ideas on the AGN ionizing spectrum and the gas den- 
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sity distribution. In particular, we also considered mixed mod- 
els with combinations of density and radiation bounded clouds 



Table 6. Metal abundances in knot 



(B96) as well as models with multiple density components. We 
spanned a very wide range of model parameters which were 
varied randomly, and selected the relatively few (about 200) 
good models which came closest to reproducing the observed 
line ratios. 

The main results are summarized in Figs. [k], [jl] which 
show the distribution of element abundances required to repro- 
duce the observed line ratios. The most remarkable feature is 
that the metal abundances are quite well constrained in spite of 
the very different assumptions made for the gas density distri- 
bution and shape of the AGN continua. In other words, models 
with very different abundances fail to match the observed lines 
ratios in knot C regardless of the AGN spectral shape and/or 
gas density distribution assumed. 

Another encouraging result is that lines from different ion- 
ization stages yield similar abundances which simply reflects 
the fact that the models reproduce the observed line ratios 
reasonably well. There are however remarkable exceptions, 
such as the [NII]/[NI] and [FeII]/[FeVII] ratios which are both 
predicted too high. Possible explanations for these differences 
have been discussed above (Sects. i^5, <L6). We stress here, 
however, that the uncertainties on [NI] have little effects 
on the derived nitrogen abundance because N + is the most 
abundant ion within the partially ionized region. Therefore, 
the best-fit N abundance decreases by only a factor ~1.3 
once the N + /N° is increased to match the observed [NII]/[NI] 
ratio (cf. Fig.JllJ). Note also that the He/H abundance is only 
poorly constrained by the models, and although models with 
He/H>0.1 are somewhat favoured, no firm conclusion about 
He overabundance can be drawn from the data. 



The derived abundances are summarized in Table |6J where 
the most striking result is the large overabundance of nitrogen 
relative to oxygen, +0.7 dex above the solar value, whose im- 
plications are discussed below. 



5.2. Comparison with other abundance estimates 

An independent estimate of metallicity can be derived from 
the measured equivalent widths of CO stellar absorption fea- 
tures, using the new metallicity scale proposed and success- 
fully a pplied to young LMC/SMC clusters by Oliva & Origlia 
(1998). In short, the method is based on the strength of the 



CO(6,3) band-head at 1.62 /im whose behaviour with metal- 
licity is modelled using synthetic spectra of red supergiants. 
The equivalent width of the stellar CO lines from the central 
100 x 100 pc 2 of Circinus are reported in Table 2 of Oliva et al. 
\\995\ and yield an average metallicity of —0.7 ±0.3, a value 
remarkably close to the oxygen abundance derived above (cf. 
Table H). 
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Element not included in the chemical evolution model 
Derived from CO stellar absorption features (cf. Sect. 5.2) 
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5.3. Nitrogen overabundance and starburst activity 

The nitrogen overabundance is of particular interest in view 
of its possible relationship with the (circum)nuclear starburst 
and N-enrichment from material processed through the CNO 
cycle. According to chemical evolutionary models of starburst 
events, the N/O relative abundance reaches a maximum value 
of [N/0]~ ±0.6 (i.e. 4 times the solar value) at about 3 10 s yr 
and remains roughly constant for several x 10 8 years (cf. Fig. 4 



of Matteucci & Padovani 1993). The nitrogen overabundance 



mostly reflects the effect of the winds from He burning red 
supergiants whose surface composition is strongly N-enriched 
by gas dredged-up from the shell where hydrogen was burned 
through the CNO cycle. The amount and temporal evolution of 
the N/O abundance depends on model details, e.g. the shape of 
the IMF and the duration of the starburst, as well as on poorly 
known parameters such as the efficiency of the dredge-up and 
the contri bution of primary N production by massive stars (e.g. 
Matteucci 1986 ). It is however encouraging to find that the ob- 
served N/O abundance (Table ^) is very close to that predicted 
at a time which is compatible with th e age of the starburst in 
Circinus (cf. Fig. 9 of Maiolino et al. 1998). It should also be 



noticed that the observed absolute abundances are about an or- 
der of magnitude lower than the model predicted values, but 
this can be readily explained if the starburst transformed only 
~ 10% of the available gas into stars, in which case the chemi- 
cal enrichment was diluted by a similar factor. This hypothesis 
is in good agreement with the fact that Circinus is a very gas 
rich galaxy (e.g. Freeman et al. 1977). 

In short, the observed N/O overabundance is fully com- 
patible with what expected for a quite old (several x 10 8 yr) 
starburst. The [NH]/Ho! and other line ratios measured in the 
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cone of Circinus are similar to those observed in many oth- 
ers Sy2's and several observational results indicate that rela- 
tively old starburst events are common in type 2 Seyferts (e.g. 
Maiolino et al. 1995). This may therefore indicate that nitro- 



gen is typically overabundant in Sy2's due to enrichment by 
the starburst associated with the AGNs. This tantalizing con- 
clusion should be however verified by detailed spectroscopic 
studies and analysis of a sufficiently large number of objects. 

6. Conclusions 

By modelling the rich spectrum of an extranuclear cloud in the 
ionization cone of the Circinus galaxy we have found that metal 
abundances are remarkably well constrained, regardless of the 
assumptions made on the shape of the ionizing continuum and 
gas distribution. This new result may open a new and inter- 
esting field of research using photoionization models to derive 
metallicities in AGNs which could in turn be related to the star 
formation activity in the recent past, i.e. old nuclear starbursts. 
In the case of Circinus, the large N/O overabundance found 
here is fully compatible with wha t exp ected from chemical evo- 
lution models of starbursts (Sect. 5.3). 

Much less encouraging are the results on the AGN ionizing 
continuum whose shape cannot be constrained by the observed 
line ratios but depends on the assumed gas density distribu- 
tion. Within the limits of the model parameters spanned here 
we somewhat favour an AGN spectrum with a "UV-bump" 
but cannot exclude that, with a different and better tuned com- 
bination of density and radiation bounded clouds, one could 
achieve similarly good fits with a power law AGN continuum 
(Sect. ph. 



We also found that photoionization models cannot repro- 
duce the observed [FeVII]/[FeII] and [NII]/[NI] ratios and ar- 
gued that these may reflect errors in the collision strengths for 
[FeVII] and rate coefficient of N-H charge exchange reactions. 
It should be noted, however, that the problem [NII]/[NI] does 
not significantly influence the derived nitrogen abundance. 

Finally, our data strongly indicate that shocks cannot play 
any important role in exciting the gas producing the observed 
line emission. 
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